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Abstract 

An extension of the Kinetic Ising model with nonuniform coupling con- 
stants on a one-dimensional lattice with boundaries is investigated, and 
the relaxation of such a system towards its equilibrium is studied. Us- 
ing a transfer matrix method, it is shown that there are cases where the 
system exhibits a dynamical phase transition. There may be two phases, 
the fast phase and the slow phase. For some region of the parameter 
space, the relaxation time is independent of the reaction rates at the 
boundaries. Changing continuously the reaction rates at the boundaries, 
however, there is a point where the relaxation times begins changing, as a 
continuous (nonconstant) function of the reaction rates at the boundaries, 
so that at this point there is a jump in the derivative of the relaxation 
time with respect to the reaction rates at the boundaries. 
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1 Introduction 



The Glauber dynamics was originally proposed to study the relaxation of the 
Ising model towards equilibrium states [lj . It is a simple non-equilibrium model 
of interacting spins with spin-flip dynamics. Kinetic models based on the Ising 
model, for example the Glauber model or the Kawasaki spin-pair exchange 
mechanism model [2J, are phenomenological models and have been extensively 
studied [THS]- It has been shown that there is a relationship between the one- 
dimensional kinetic Ising model at zero temperature and diffusion annihilation 
in one dimension [3 ]. In [4 , using a damage spreading method, the sensitivity of 
the time evolution of a kinetic Ising model with Glauber dynamics against the 
initial conditions has been investigated. The full time dependence of the space- 
dependent magnetization and of the equal time spin-spin correlation functions 
were studied in [5] . Non-equilibrium two-time correlation and response functions 
for the ferromagnetic Ising chain with Glauber dynamics have been studied 
in [HE]. 

Combinations of the Glauber and the Kawasaki dynamics have been also 
considered [8HTU]. Most studies are focused on uniform lattices where reaction 
rates are site-independent. It is known that the ordinary Glauber model on a 
one-dimensional lattice with boundaries at any temperature, shows a dynamical 
phase transition [11]. The dynamical phase transition is controlled by the rate 
of spin flip at the boundaries, and is a discontinuous change of the derivative 
of the relaxation time towards the stationary configuration. Among the sim- 
plest generalizations beyond a uniform system is a lattice with alternating rates. 
In [12Ml4| . the steady state configurational probabilities of an Ising spin chain 
driven out of equilibrium by a coupling to two heat baths have been investi- 
gated. An example is a one-dimensional Ising model on a ring, in which the 
evolution is according to a generalization of Glauber rates, such that spins at 
even (odd) lattice sites experience a temperature T e (T ). In this model the 
detailed balance is violated. The response function to an infinitesimal magnetic 
field for the Ising-Glauber model with arbitrary exchange couplings has been 
studied in [15] . Other generalizations of the Glauber model consist of, for exam- 
ple, alternating- isotopic chains and alternating-bound chains ( |16j for example). 
In |17j . an asymmetric generalization of the zero-temperature Glauber model 
on a lattice with boundaries was introduced. There it was shown that in the 
thermodynamic limit (when the lattice becomes infinite) the system shows two 
kinds of phase transitions. One of these is a static phase transition, the other 
a dynamic one. The static phase transition is controlled by the bulk reaction 
rates, and is a discontinuous change of the behavior of the derivative of the 
stationary magnetization at the end points, with respect to the reaction rates. 
The dynamic phase transition is controlled by the spin flip rates of the particles 
at the end points, and is a discontinuous change of the relaxation time towards 
the stationary configuration. Other phase transitions induced by boundary con- 
ditions have also been studied ( [I8l - l20] for example). Another generalization of 
the Glauber model was introduced in [3T]. In this generalization, the processes 
are the same as those of the ordinary Glauber model, but the rates depend on 
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three free parameters, rather than one free parameter as in the ordinary Glauber 
model. Finally, this model was further generalized to the case where the number 
of interacting sites is more than three and the number of states at each site is 
more than two. This model too violates detailed balance. In [35], an inherent 
spin anisotropy (kinetic disorder) in the Glauber- Ising model was introduced. 
It was shown both analytically and numerically, that there is a slow logarithmic 
factor in the decay of the density of kinks at large times. In [23], the problem of 
the effect of quenched impurities on one-dimensional non-equilibrium Glauber- 
Ising-type models has been investigated numerically. It was shown that the 
model has a continuous phase transition to an absorbing state. Also a mean- 
field approach has been used to study the Glauber-type stochastic dynamics 
of a model on a square lattice in which two interpenetrating square sublattices 
have spins that can take two values, alternated with spins that can take the four 
values [2"3] . 

The behavior of an Ising model with nonuniform coupling constants at low 
temperatures on a one dimensional periodic lattice has been discussed in |25j . 
The static behavior of an Ising model with nonuniform coupling constants on 
a one-dimensional lattice with boundaries was investigated in |26j . Detailed 
balance was used in [2"51l2l)] to propose reaction rates for the system. In [53] the 
time- independent solution to the evolution equation of the expectation values 
of spins was studied. This solution was expressed in terms of a transfer matrix. 
While it is true that the ordinary Ising model does not exhibit any phase tran- 
sition in finite temperatures (the expectation values of the spins vanish if there 
is no external magnetic field), this is not necessarily the case for the Ising model 
with inhomogeneous boundary conditions. It was shown that in the thermody- 
namic limit different phases could occur for this system, according to whether 
the eigenvalues of the transfer matrix are less than or larger than one. 

In this paper the dynamical properties of kinetic Ising model with nonuni- 
form coupling constants are investigated. The nonuniformity could arise from 
either a nonuniformity in the coupling constants or a nonuniformity in the tem- 
perature. However, as long as the dynamics of the system is considered, only the 
ratios of the coupling constants and the temperature are important, so without 
loss of generality one can assume the temperature to be constant and put all the 
nonuniformity in the coupling constant. So a nonuniformity in the (effective) 
coupling constants could be due to a nonuniformity in the temperature. It could 
also be due to a nonuniformity in the interaction between neighboring spins, or 
due to a nonuniformity in some sort of inertia, which opposes the changes in 
the state of the system. 

It is shown that there are cases where the system exhibits a dynamical 
phase transition. There are two phases: the fast phase and the slow phase. For 
some region of the parameter space, the relaxation time is independent of the 
reaction rates at the boundaries. This is the fast phase. Changing continuously 
the reaction rates at the boundaries, there is a point where the relaxation time 
begins increasing, so that at this point there is a jump in the derivative of the 
relaxation time with respect to the reaction rates at boundaries. This is the 
dynamical phase transition and the region where the relaxation time changes 
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with reaction rates at boundaries is the slow phase. So the dynamical phase 
transition studied here, is a discontinuity in the derivative of the relaxation time 
(from zero to nonzero) with respect to reaction rates at the boundaries. 

The scheme of the paper is as follows. In section 2, the model is introduced, 
and the evolution equation for the spin expectation values is obtained. In section 
3 the relaxation of the spin expectation values towards their stationary state is 
investigated, through a study of the time scales of this relaxation, and a general 
equation is obtained for these time scales. In section 4, some examples are 
studied in more detail. It is seen that some of these show two phases, the fast 
phase and the slow phase. Section 5 is devoted to the concluding remarks, and 
section 6 is an appendix on the site-link notation and the relaxation time. 



2 One-dimensional Ising model with nonuniform 
coupling constants 

Consider a one-dimensional lattice with (L + 1) sites, labeled from to L. At 
each site i, there is a spin variable Sj, which could be +1 (spin up, f), or —1 
(spin down, \). These spins interact through the Hamiltonian, 



H = H 



where 1-L a is the Hamiltonian interaction for the link a: 

J a is the coupling constant in the link a, and 

,:=\- 0) 

The link a links the sites a — fi and a + fi, so that a ± fi are integers, and a 
runs from fi up to (L — fi). Throughout this paper, sites are denoted by Latin 
letters which represent integers, while links are denoted by Greek letters which 
represent integers plus one half (ji), so that the link a joins the sites (a — fi) 
and (a + fi), while the site i joins the links (i — fi) and (i + fi). TL' and W L 
correspond to interactions at the boundaries. A more extended explanation of 
the site-link notation is presented in the appendix. 

A nonuniform Glauber model gives the dynamics of the Ising model with 
nonuniform coupling constants ( J a ) such that the rate of a spin flip is determined 
through its interaction with its two neighboring sites and a heat bath at the 
temperature T. Introducing 

K - : = & < 4 > 

where &b is the Boltzmann's constant, it was shown in [2 6) that assuming nearest 
neighbor interaction, and that in each step only one spin flips, detailed balance 
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gives 

u[(S', 8j ) -> (5', - Si )] = [1 - « 3 - tanh(^_ M Sj _! + s j+1 )], (5) 

where S" denotes the configuration of the lattice apart from site j, and Tj's 
are constants. So that the spin at the site j flips according to the following 
interactions with the indicated rates. 

ttt^tlt and 4- 4- 4- 4- 1 4- with rate 1 - tanh(X i _ AI + K j+fl ), 
tlt^ttt and 1 1 4- ->■ I 4- 4- with rate 1 + tanh(X i _ Al + Kj +fl ), 
1 1 4 1 4- 4- and llt^ltt with rate 1 - tanh^-^ - Kj +fl ), 
It t~> lit and tll^ttl with rate 1 + tanh^-^ - K j+fll ), 

(6) 

where Tj's have been taken independent of j, and set to one by rescaling the 
time, as in [25j 26 . At the boundaries, there are other interactions as well 
(corresponding to TL' Q and 7i' L ). The spin of the zeroth site may flip like 

1 1 — y 4- 4- with rate gi , 

tt^tt with rate g 2 , 

It^tt with rate g 3 , 

I I ->■ 1 4 with rate g 4 , (7) 
and the spin of the L'th site may flip like 

It - ^4-4- with rate hi, 

tt^tl with rate hi, 

t4-->-tt with rate h 3 , 

| 4- -> 4- 1 with rate /14. (8) 

It is known that with the rates ([6]) , the time derivatives of the one-point functions 
in the bulk are expressed in terms of only the one-point functions. To make this 
true for the boundaries as well, the following relations should hold [TT] . 

9i +9i = 92 +93, 

hi + hi = h 2 + h 3 . (9) 

The first relation, for example, means that the sum of conditional spin flip rates 
at the zeroth site is independent of the state of the first site. This does not, 
however, means that the flip rate at the zeroth site is independent of the state 
of the first site; since, for example, (gi + g 4 ) itself is not a (conditional) spin flip 
rate. 

Thus the evolution equation for the expectation values of the spins are 

(sj) = - 2 ( Sj ) + [tanh(Xj-_ p + K j+fl ) + tanh(if J _ M - K j+tl )} 

+ [tanh( J ftT i _ AI + Kj +tl ) - tanh^-^ - K j+I1 )] (s j+ i), < j < L 
(so) = (.93 - 9i) - (32 + 53) (so) + (51 - 92) (si), 

(sl) = {h 3 - hi) - (h 2 + h 3 ) (s L ) + {hi - h 2 ) (s L -i). (10) 
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The static solution ((s) s t) was studied in [33]. Here the dynamics (the relaxation 
towards the static solution) is addressed. 



3 The relaxation of the system towards equilib- 
rium 



The homogeneous part of eq. ([TO]) can be written as 

fr) = h) ( Sl ). 



(11) 



A brief introduction regarding the relaxation time of the system has been con- 
tained in the appendix. Denoting an eigenvalue of h by E, and the corresponding 
eigenvector by x, it is seen that 



Exj — —2xj + [tanh(i£j_ M + K ]+ ^) + tanh(i4T J _ Al - K ]+ ^)] Xj-i 

+ [taxJa.(Kj^ fl + Kj + n) - tanh(K J _ A( - Kj +fl )} x J+u < j < L 



Ex Q = - (g 2 + 93) x + (gi - g 2 ) x%, 
Ex L = - (h 2 + h 3 ) x L + (hi - h 2 ) xl-i, 



which can be written as 
where 



X a :— 



u a—fj. 



> Xa+f, 



and 



D: 



( 







1 



sinh(2 Kj-fj) cosh(2 Kj-^) + cosh(2 K 



1 ~~\~~ — 



\ 

E 

2JJ 



V wah.(2K j+h ) sinh(2JC J+Al ) 
Using the recursion relation (fT5|i repeatedly, one arrives at 

X a — D a p Xp , 

where 

D a /3 '■— Da-n ■ ■ ■ Dp +ll 
(fT3"|) and (|14l) are boundary conditions for X, and can be rewritten as 




(12) 
(13) 
(14) 

(15) 
(16) 



(17) 



(18) 
(19) 

(20) 
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and 



oi- 



and 



^i_A (h 2 + h 3 + E\ 

|. (21) 

x L V h x - hi 



(22) 




Defining 

it is seen that 



(23) 



W := (fta - hi h 2 + h 3 + E), (24) 



WI £ -, = 0. (25) 
Then, using (IT51) for a = L — /j, and /3 = /j,, one arrives at 

W^l-^ = 0, (26) 

where 

/ .9i ~ 92 \ 
V:=[ . (27) 

\92+93, + E) 

(|2^|) is a polynomial equation of order (L + 1) for E, which has (L + 1) roots 
for E. All roots should have nonpositive real parts. The root with largest real 
part determines the relaxation time of the system. The aim is to investigate the 
behavior of this root in the thermodynamic limit (L — > oo). 

4 Special cases 

1 Uniform coupling constant: 

K a - K, (28) 

This is the same kinetic Ising model, the dynamical phase transition of 
which induced by boundary terms was investigated in |17j , using a different 
method. To fix notation and become familiar with the transfer matrix 
method, let's use this method for this case. In this case Dj is independent 
of j: 

D :=\ ' ] . (29) 

\-l (2 + E) coth(2K)J 

Denoting an eigenvalue of the the matrix D by z, it is seen that 

E = -2+(z + Z - 1 ) tanh(2iT). (30) 



G 



So the product of the eigenvalues is one. (flZB")) recasts to 

WD L ~ 1 V = 0. (31) 



D can be written as 



I 1 




D = -^~. I II II I- 



from which one obtains 
1 



^ = 7^- | • (33) 

« \ L-l_ z -i+l _ z i + z -i 




z L-2_ z -L+2 _ Z L-1 + Z -L+^ 



V = 0. (34) 



Substituting E, using (f3T)|). one arrives from the above equation at a poly- 
nomial equation of order (2 L + 4) for z. Two roots of this equation are 
±1, which are not the roots of (|3"Tj) . as (f3~4"|) has been obtained from (|31l) 
by multiplying it by (z _1 — z). The other (2L + 2) roots consist of (L + 1) 
pairs. In each pair the two roots are inverses of each other, and both roots 
of each pair give the same value for E. 

The roots of (|34[) for z are either unimodular or not. For a unimodular 
solution 

z=:e ie , (35) 

equation (fM)) becomes 



/sin(£ -2)9 - sin(L - 1) 
W | 7 = 0, (36) 

\sin(£-l)0 -sinL0 

which can be written like 

F(B) = 0. (37) 

It is seen that in the thermodynamic limit, the entries of the matrix ap- 
pearing in the left-hand side of (|3"6"]) change sign when the value of 9 is 
changed by (ir/L), while V and W do not change. So the sign of F{ff) is 
opposite of the sign of F[9+ (ir/L)], meaning that for any value of #o there 
is at least one root for 9 between 9q and [9o + (ir/L)]. So, in the thermo- 
dynamic limit all of unimodular z's are roots of |34|) . These correspond 
to real values of E, with the maximum 

E a = -2 + 2 tanh(2 1^1), (38) 
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which corresponds to the relaxation time 

1 

T ° ~ 2-2tanh(2|if|)' (39) 

If there are no nonunimodular solutions for z, this is the relaxation time of 
the system, which does not depend on the reaction rates at the boundaries. 
Let's call this the fast phase. If, however, there is a value for E the real part 
of which is larger than the right-hand side of (|3"51) , then the relaxation time 
is larger than tq and, as it will be seen, does depend on the boundaries. 
This is called the slow phase. 

If there is a nonunimodualr root for z, then there is a root z with modulus 
larger than one. Then, in the thermodynamic limit (|34[) becomes (for this 
root) 

-1\ 

W I I F = 0, (40) 

or 

[{h 2 - h^z- 1 + h 2 + h 3 +E] [(g 2 - gi ) z' 1 + g 2 + g 3 + E] = 0. (41) 

Using (|3T)]) . this equation turns out to be 

{[h 2 - hi + tanh(2 iT)] z^ 1 + [tanh(2iT)] z + h 2 + h 3 - 2} 
x{[52 ~ 9i + tanh(2 K)] z" 1 + [tanh(2 K)] z + g 2 + g 3 - 2} = 0, (42) 

which is exactly the same as equation (27) in [T7]. So the rest of discussion 
is exactly similar to that of [17] . The root obtained from (l42t does depend 
on the reaction rates at the boundaries. Such a root corresponds to a 
relaxation time larger than To, iff 

[Reiz + z' 1 )} [sgn(*0] > 2. (43) 

This gives the following coexistence surface for the two phases 

52 + 53 + [sgn(X)] (32 - 3i) - 2 [1 - tanh(2 \K\)] = 0, (44) 

or 

[4 tanh(2 \K\) + 2 - g 2 - g 3 ] [tanh(2 \K\) - gi + g 2 ] 
+ [tanh(2|X|)](3 2 + 33 - 2) = 0, (45) 

and similar solutions with hi's substituting g^s. 

In summary, for some values of reaction rates at the boundaries the relax- 
ation time is To , which is independent of the reaction rates at the bound- 
aries. At a certain point, however, the relaxation time starts increasing 
from tq and being dependent on the reaction rates at the boundaries. The 
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dynamical phase transition is this discontinuous change of the derivative 
of the relaxation time towards the stationary configuration. 

As examples for this transition, let's consider two cases. Regarding the 
coexistence curve (T4"4")l . it is seen that for 



91 

94 



52, 
53, 



one has 




53 < 2[1 - tanh(2 \K\ 
53 > 2[1 - tanh(2 \K\ 



slow phase 
fast phase 



(46) 
(47) 



Regarding the coexistence curve (|4"5|) . it is seen that for 



.9i = .93, 

94 = 52, 

.92 + .93 = 2, 



one has 



.93 — 52 < tanh(2 \K\), slow phase 
.93 — 52 > tanh(2 \K\), fast phase 



(48) 
(49) 



2 K n 



K', a < L' 
K, a > V 



It has been assumed that in the thermodynamic limit both L' and [L — L') 
tend to infinity. Also, without loss of generality it can be assumed that 
\K\ > \K'\. This is a lattice consisting of two parts, on each part the 
coupling constant is uniform. Equation (|26[) becomes 



(50) 



W f)L-L -1 jj,L -1 y _ q 



where D is defined through (|29|) . D' is similar to D but with K replaced 
by K', and 



Drj : = 







1 



sinh(2 K') cosh(2 K') + cosh(2 iiT) 



E 

1+ 2 



(51) 



\ sinh(2 K) sinh(2 X) 

Denoting an eigenvalue of D (D') by z (z 1 ), it is seen that holds and 

£ = -2 + (z' + z'" 1 ) tanh(2X')- (52) 

Equations (|50|) . (f30|) . and (l52t are used to obtain the values of z and z'. 
First assume that z is unimodular. Then E is real and from (|52[) it is 
seen that z' is either unimodular or real. If z' is unimodular too, then the 
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possible values for z are not all of the points of the unit circle. Defining 9 
by (|35|) . it is seen that those values of 9 which correspond to unimodular 
values for z' are 



cos 1 



tanh(2|Jf'|) 
tanh(2|if|) 



< \9\ < 7T — COS 



tanh(2|iT'| 



tanh(2 \K\) 



(53) 



For each value of 9 out of this region, there are two real values for z' which 
are inverses of each other. For that value of z' the modulus of which is 
greater than one, one arrives from (1501) at 

/ sin(X — L' — 2)0 — sm(L — L' — 1) 6\ iz'- 1 -l\ 

W \ M L ' ) v = o, 

\sin(L-iy- 1)6* -sm(L-L')6 J \ 1 -z'J 

again showing that if (L — L') tends to infinity, all of the (remaining) 
values of 9 are roots of this equation. So one piece of the solutions for E 
corresponds to the full unit circle for z and the full unit circle plus two real 
segments for z' . The relaxation time corresponding to these eigenvalues is 
To from (|3"5|) . It is seen that it is the larger coupling which determines the 
relaxation time (provided the part of the lattice which corresponds to this 
coupling becomes infinite). If all of the solutions for z are unimodular, tq 
is the relaxation time of the system and the system is in the fast phase. 

The system would be in the slow phase (where the relaxation time is larger 
than To and does depend on the rates at the boundaries) iff there exists a 
solution for z for which P3"| holds. In this case, one arrives at 

w ( 1 b v i l j v ^°- (55) 

where z and z' have are solutions which have moduli larger than one. One 
then arrives at 

{[h 2 -hx + tanh(2 K)] z^ 1 + [tanh(2 K)] z + h 2 + h 3 - 2} 
x[z sinh(2iT) + z sinh(2 K')\ 

x {[92 - 91 + tanh(2 K')] z'- 1 + [tanh(2 K')] z' + g 2 + 93 - 2} = 0. (56) 

From (|3T)|) and (|52p . it is seen that the real parts of [z tanh(2if)] and 
\z' tanh(2 K')], and hence the real parts of [z sinh(2 A')] and \z' sinh(2 K')] 
have the same sign, or both are zero. So the second factor in the left hand 
side of (|56[) does not vanish unless z and z' are (±i). But these values 
of z and z' arc unimodular hence not acceptable, as the roots of ([56"]) . 
So in order that ([56]) is satisfied, either the first or the third factor in 
the left hand side of (f56|) should vanish. This is similar to (|42|) . with the 
couplings corresponding to each part of the lattice instead of the uniform 
coupling. One then arrives at (|44|) and (|45|) with K replaced by K', or 
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similar expressions with K' and t^'s replaced by K and hi's respectively, 
as the coexistence surface. Of course from the possible roots of (|56"T) . that 
root is chosen which corresponds to the larger relaxation time (or larger 
real part for E) , and one should also check that this root satisfies fH]) . 

To summarize, in the fast phase the relaxation time is determined by 
the coupling with larger modulus. In the slow phase the two ends of the 
lattice behave independently, the relaxation time being determined by the 
largest. 



3 K a = 



K>, 
K" 



a<L' 

L' < a < L - 

L - L" < a 



L" 



It has been assumed that in the thermodynamic limit both L' and L" 
tend to infinity, while (L — L' — L") remains finite. This is a lattice with a 
finite part inside which, where outside that part the coupling is uniform. 
Equation ([26"]) becomes 



wb" L "- 1 D L . 



L"+n L'-fj, 



f)'L'-l y = q 



(57) 



where D' and D" are defined through (|29|). but with K replaced by K' 
and K" respectively. Similar to the previous example, one defines z' and 
z" similar to (1301) and (|52[) . Again similar to the previous example, there 
is a relaxation time corresponding to the fast phase, where from z' and z" 
at least one is unimodular, which is obtained from (|39|) . with \K\ being 
the larger between \K'\ and \K"\. In order that the system be in the slow 
phase, there should be a solution where both z' and z" are nonunimodular. 
As for any value of E, the product of the eigenvalues of bi is one, one 
can take z' and z" so that both have moduli larger than one. One then 
arrives at 

W 



-'-i 



D, 



' L — L" -\-fj, L' — jj, 

J>1 x i 



V = 0, (58) 



or 



(1 -•:") Dl l- +i ,l> „ j j [(-z'- 1 1) V] =0. 

(59) 

Putting the first or the third factor equal to zero, a result similar to that 
of the previous example is obtained. But there may be another solution 
which is obtained by letting the second factor vanish: 











w 

















1 -z") D 



h—h" L' — jj, 



= 0. 



(60) 
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This new solution (if it exists) can change the relaxation time. If the 
relaxation time corresponding to the solution to (|60| is larger, then this 
relaxation time is the relaxation time in the slow phase. So, the introduc- 
tion a finite part in the lattice can change the relaxation time and hence 
the coexistence surface. In fact, qualitatively this finite part plays the role 
of new boundaries introduced in the lattice. 

4 Alternating Coupling constant K a = (— l) a ^^ K 
Defining 



'0 



D r 



-(2 + E) coth(2 iT) y 
1 

1 (2 + E) coth(2 K) . 



it is seen that (f2"S|) becomes 

W Dl- r (D e A>) (L ~ 2+r)/2 V = 0, 
where r is (1) when L is even (odd). Defining z through 



E 



-2+ (z + z^ 1 ) tanh(2iT), 



it is then seen that 



1 -(z + z- 1 ) 

K {z + z- x ) I -(z + z- 1 ) 2 
showing that the eigenvalues of (D c D ) are (— z 2 ) and (— z~ 2 





(61) 

(62) 
(63) 

(64) 



So, 



(D e D y = 



("I) 9 



1 





V 



(65) 



1 



Again (|6^|) shows that tq is the relaxation time of the system (the system 
is in the fast phase), if all of the values of z are unimodular. To have a 
larger relaxation time (the system be in the slow phase) , there should be 
a value of z with modulus greater than 1. In fact, there should be a value 
of z satisfying ()43[) . In that case, in the thermodynamic limit (q —5- oo) 
one arrives from (|62p to 







WDl~ r 


)i 







-i) v]=o, 



(66) 
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or 



'(-1) 



W 



l-r 



1 



-1) V] = 0. 



(67) 



It is seen that for r = 1 (odd L), this is exactly ([4"2j). For r = (even L), 
the second factor in the left-hand side of ([66]) is the same as the second 
factor in (1421) . while the first factor is the first factor in (|42|) with 2 and 
AT replaced by (—z) and (—if), respectively. So the coexistence surfaces 
(|44p and (|45p are recovered here as well, while the coexistence surfaces 
corresponding to h^s are the same as those obtained for the uniform lattice 
if L is odd, and the same with A replaced by (—A) if L is even. 

One could also arrive at this result by making a correspondence between 
this case (alternating coupling constants) and case 1 (uniform coupling- 
constants). The analogy goes as the following. 



K a 

(91, h, 9a, fa) 
(hi,h 2 , ha, hi) 



[j/2] 



(-1) 

(-1) C 
(51,52,53,54), 

J (hi, h 2 , ha, hi) 
[ (h-2,hx,hi,ha) 



for odd L 
for even L 



(68) 



where [y] means the largest integer not greater than y. It is seen that 
in this way, if the untilded varibles and constants satisfy (|10p till (|T2j) . 
the tilded varibles and constants satisfy (fT0|) till (fi"2|) as well, and if A Q is 
alternating, K a is obviously uniform. The meaning of this correspondence 
between Xj's and Xj's is that, for example, xq and x\ are kept fixed, while 
X2 and X3 have been changed to (—x 2 ) and (—Xa), respectively, and A]/ 2 
and A'5/2 have been kept fixed, while -K3/2 has been changed to (—^3/2)- 
Obviously the expressions (ifi/2 ^0 ^1), (K 3/2 x\ x-i), (^5/2^2^3) have 
all remained intact, while the coupling is now uniform, if it has been 
alternating. If L is odd, no change in the boundary rates is required and 
the coexistence curves would be exactly the same as those found in the case 
of uniform coupling. If L is even, that part of the coexistence curve which 
arises from the g's is still intact, as <?'s have remained intact. But the part 
due to /i's is changed. And the change comes from the fact that (h 2 — h%) 
has been changed to (h± — h 2 ), while (h 2 + ft-3) has remained intact. One 
can see that in the equations of the coexistence curves involving h's ((|44p 
and (|45|) but with g's replaced with h's), this change is equivalent to a 
change of z to (—z) and K to (-K). 



5 Concluding remarks 

A one dimensional Ising model with nonuniform coupling constants on a lattice 
with boundaries was studied. Detailed balance was used to obtain the evolu- 
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tion of this model. The relaxation of the resulted system towards its stationary 
solution was studied. A general formulation was obtained to determine the 
relaxation times and possible dynamical phase transitions, where the system 
moves from the fast phase (with the relaxation time independent of boundaries) 
to the slow phase (with the relaxation time depending on boundaries). Some 
special examples were studied in more detail, including cases where the lattice 
consists of two infinite homogeneous parts and one middle (possibly nonhomo- 
geneous part), where it was shown that this middle part can induce further 
dynamical phase transitions, behaving like an effective boundary. 

Part of the results obtained here is similar to those corresponding to the 
case of uniform coupling constants. The similarity is that the system is in the 
fast phase when the boundary couplings are high enough so that it is the bulk 
reaction rates that determine the relaxation, and goes to the slow phase when 
the boundary reactions are less than some critical value. Another similarity 
is that each boundary behaves essentially independent to the other boundary. 
When the system consists of large parts which are essentially uniform, the new 
feature in the relaxation is that the system behaves as though it is consists 
of independent systems each having two boundaries. So for such systems new 
coexistence curves arise. There are of course cases where no large uniform blocks 
are present. Then the only similarity with the uniform case is that increasing 
the boundary rates makes them unimportant in the relaxation time. 

6 Appendix 

6.1 notation for a lattice 

A directional lattice is a collection of sites and directional links, with a relation 
between links and sites. A link 1 has a negative boundary (denoted by d~l), and 
a positive boundary (denoted by d + l). The relation of the site i to the link 1 is 
that either i is d~l, or i is or i is not in 1. To each site i there corresponds 
a set of outgoing links S + (i), and a set of incoming links S~(i): 

S+(i) :={l|0-l = i}, 

S~(i) := {1 | 8+l = i}. (69) 

In a one dimensional (connected) lattice these relations become simple. The 
lattice sites can be denoted by integers from a subset Q of the integers, with 
the property that if an integer is between two members of Q, it is a member of 
Q itself. Q can be bounded from either above or below. If it is bounded from 
above, then the lattice has a positive boundary, denoted by d + Q, which is the 
largest member of Q. Similarly the smallest member of Q (if it exists) is the 
negative boundary of Q which is denoted by d~Q. S + (i) consists of exactly one 
link, unless i is d + Q, in which case S + (i) is empty. Similarly, S~(i) consists 
of exactly one link, unless i is d~Q, in which case S~(i) is empty. So one 
can denote the outgoing and incoming links corresponding to the site i, simply 



14 



through 

i±fj,:=S ± (i), (70) 

from which one also has 

d ± a = a±fi. (71) 

If the lattice is closed and has L sites, one can still use the above notation, 
provided an equivalence between b and (b + L) is assumed, where b can be a site 
(an integer) or a link (an integer plus half). 

6.2 the relaxation time 

Denoting the deviation from the static configuration of the system by x, it is 
seen that x satisfies 

y = hy, (72) 

which is nothing but the closed form of (fTTj) . To solve this, one expresses y as 
a linear combination of the (generalized) eigenvectors of h: 

y(t) = ^2c x (t)xx, (73) 

A 

where x\ is the (generalized) eigenvector of h corresponding to the eigenvalue 
E\. If all of the generalized eigenvectors are eigenvectors, evolution equations 
for c\s decouple as 

ca=#aca, (74) 



resulting in 

So one arrives at 



c x (t) =c A (0) cxp(E x t). (75) 
y(t) =J2 ca(0) eMExt)x x . (76) 



Initial condition determines ca(0)'s. It is seen that for large times the leading 
term of the right hand side of (l76l) is that term which corresponds to the largest 
real part of Ex, unless the initial condition is fine tuned so that ca(0) vanishes 
for that A. So the relaxation time is 

' (77) 



max[Re(i?A)] ' 

In the thermodynamic limit, it may happen that the maximum of Re(_EoO tends 
to zero, so that the relaxation time goes to infinity, meaning that the relaxation 
is no longer exponential but say power law. But it may also happen that this is 
not the case, and the maximum of Re(E\) tends to a negative number, in which 
case the relaxation time remains finite even in the thermodynamic limit. 
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